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We calculate QCD corrections to transversely polarized Drell-Yan process at a mea- 
sured Qt of the produced lepton pair in the dimensional regularization scheme. The Qt 
distribution is discussed resumming soft gluon effects relevant for small Qt- 
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Hard processes with polarized nucleoli beams enable us to study spin-dependent dy- 
namics of QCD and the spin structure of nucleon. The helicity distribution Aq(x) 
of quarks within nucleon has been measured in polarized DIS experiments, and 
AG(x) of gluons has also been estimated from the scaling violations of them. On 
the other hand, the transversity distribution 5q(x), i.e. the distribution of trans- 
versely polarized quarks inside transversely polarized nucleon, can not be measured 
in inclusive DIS due to its chiral-odd nature™ and remains as the last unknown 
distribution at the leading twist. Transversely polarized Drell-Yan (tDY) process is 
one of the processes where the transversity distribution can be measured, and has 
been undertaken at RHIC-Spin experiment. 

We compute the 1-loop QCD corrections to tDY at a measured Qt and az- 
imuthal angle <f> of the produced lepton in the dimensional regularization scheme. 
For this purpose, the phase space integration in D-dimension, separating out the rel- 
evant transverse degrees of freedom, is required to extract the cxcos(20) part of the 
cross section characteristic of the spin asymmetry of tDYPThe calculation is rather 
cumbersome compared with the corresponding calculation in unpolarized and lon- 
gitudinally polarized cases, and has not been performed so far. We obtain the NLO 
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0(a s ) corrections to the tDY cross section in the MS scheme. We also include soft 
gluon effects by all-order resummation of logarithmically enhanced contributions 
at small Qt ("edge regions of the phase space") up to next-to-leading logarithmic 
(NLL) accuracy, and obtain the first complete result of the Qt distribution for all 
regions of Qt at NLL level. 

We first consider the NLO 0(a s ) corrections to tDY: h\(P\, si) + h<z(P2, S2) — > 
l(k\)+l(k2)+X , where hi, hi denote nucleons with momentum Pi, P2 and transverse 
spin Si,s<2, and Q = ki + ki is the 4- momentum of DY pair. The spin dependent 
cross section A^dcr = (da(si, S2) — da(si, — S2))/2 is given as a convolution 

At da = J dxidx 2 6H(xi , X2] ^f) ^rda(si : S2; ^f), (1) 
where fiF is the factorization scale, and 

5H(xi , x 2 ;hf) = ^2^[Sqi(xi; /j, F )Sqi(x 2 ; /j,f) + Sqifai; fXF)8qi(x 2 ; Hf)] (2) 

i 

is the product of transversity distributions of the two nucleons, and A-r<i<T is the 
corresponding partonic cross section. Note that, at the leading twist level, the gluon 
does not contribute to the transversely polarized process due to its chiral odd nature. 
We compute the one-loop corrections to A^dcr, which involve the virtual gluon 
corrections and the real gluon emission contributions, e.g., q(pi,si) + q(p2,S2) — * 
l(ki) + l(k2)+g, with pi = XiPi. We regularize the infrared divergence in D = 4 — 2e 
dimension, and employ naive anticommuting 75 which is a us ual p rescription in the 
transverse spin channelPIln the MS scheme, we eventually get to NLO accuracy, 

A ' l] " N cos (20) [X (Q T , Q 2 , y) + Y (Q T , Q 2 , y)] , (3) 



dQ 2 dQ T dyd(f) 

where N — a 2 / (3 N c S Q 2 ) with S = (Pi + P2) 2 , y is the rapidity of virtual photon, 
and 4> is the azimuthal angle of one of the leptons with respect to the initial spin 
axis. For later convenience, we have decomposed the cross section into the two parts: 
the function X contains all terms that are singular as Qt — > 0, while Y is of 0(a s ) 
and finite at Qt = 0. Writing X = + X^ as the sum of the LO and NLO 
contributions, we havJ33 1(0 = SH(x° 1 , x° 2 ■ up) S(Q 2 t ), and 



xM = ^-Cf\6H( X 1 X °;»f) 2(^p^) - + (_ 8 + **) 



1 




+ S(Q T )hx 



(Q T )+ Mf 



) + (4 - xD 



'(4) 



where x\ — \fr e y , x\ = \pr e~ v are the relevant scaling variables with r = Q 2 /S, 
and SPqq (z) = 2z/(l — z) + + (3/2) 6(1 — z) is the LO transverse splitting function.^ 
In 0, the terms involving 5(Q 2 ^) come from the virtual gluon corrections, while 
the other terms represent the recoil effects due to the real gluon emissions. For the 
analytic expression of Y , see RefH Eq. J2J gives the first NLO result in the MS 
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scheme. We note that there has been a similar NLO calculation of tDY cross section 
in massive gluon scheme.® We also note that, integrating J3J over Qt, our result 
coincides with the corresponding Qr-integrated cross sections obtained in previous 



works employing massive gluon scheme2land dimensional reduction scheme ® via 
the scheme transformation relation^ 

The cross section (jSJ becomes very large when Qt <C Q, due to the terms be- 
having ~ a s ln(Q 2 /Qt)/Qt an< ^ ~ &s/Qt m the singular part X. It is well-known 
that, in unpolarized and longitudinally polarized DY, large "recoil logs" of simi- 
lar nature appear in each order of perturbation theory as a™ In 2 ^ 1 (Q 2 /Q T )/Q T , 
a™ ln 2n ~ 2 (Q 2 /Q T )/Q T , and so on, corresponding to LL, NLL, and higher contri- 
butions, respectively, and that the resummation of those "double logarithms" to all 
orders is necessary to obtain a well-defined, finite prediction of the cross section.® 
Because the LL and NLL contributions are universalpl we can work out the all- 
order resummation of the corresponding logarithmically enhanced contributions in 
© up to the NLL accuracy, based on the general formulation^ of the Qt resum- 
mation. This can be conveniently carried out in the impact parameter b space, 
conjugate to the Qt space. As a result, the singular part X of |J5J is modified into 
the corresponding resummed part, which is expressed as the Fourier transform, 

X E e i db h -MbQ T )e s ^\C qq <g> %) U, j|) (C qq ® Sqt) (x», |) 

+{xl <- (5) 

Here b = 2e fB , and the large logarithmic corrections are resummed into 
the Sudakov factor e s ^ with S(b,Q) = - J^ /b2 (dK 2 / K 2 ){(ln ^)A q (a 8 ( K )) + 
-B g (a a («))}. The functions A q , B q as well as the coefficient functions C qq ,C qq 
are cal culable in perturbation theory, and at the present accuracy of NLL, we 
get!5EU,(a fl ) = (a s /7r)C F + (a s /2^) 2 2C7 F {(67/18- 7r 2 /6)C G - 5iV//9}, B q {a s ) = 
-3C F (a s /2Tr), C qq (z,q t )= C qq (z,a s ) = 5(1 - z){l + (a s /4n)C F (7r 2 - 8)}. We 



have utilized a relation" 1 between A q and the DGLAP kernels in order to obtain 
the two- loop term of A q . The other contributions have been determined so that 
the expansion of the above formula 10) in powers of a s (/ii?) reproduces X of Qt, 
(@J to the NLO accuracy. Eq. J3J with (JJJ presents the first result of the NLL Qt 
resummation formula for tDY. The NLO parton distributions in the MS scheme 
have to be used. 

One more step is necessary to make the QCD prediction of tDY. Similarly to 
other all-order resuumation formula, our result (JjjJ is suffered from the IR renor- 
malons due to the Landau pole at b = (6 /Q)e ( - 1 / 2/3oQs ^' ) - ) in the Sudakov factor, 
and it is necessary to specify a prescription to avoid this singularity. Here we de- 
form the integration contour in (jSJ i n t he complex b space, following the method 

in t\i 

introduced in the joint resummation. Obviously prescription to define the b in- 
tegration is not unique reflecting IR renormalon ambiguity, e.g., "6* prescription" 
to "freeze" effectively the b integration along the real axis is frequently used!^ The 
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Fig. 1. Q T distribution A T da/dQ 2 dQ T dyd<f> at = 100 GeV, Q = 10 GeV, y = 4> = 0. 

renormalon ambiguity should be eventually compensated in the physical quantity 
by the power corrections ~ (6Aqcd)™ (n — 2,3,...) due to non-perturbative ef- 
fects. Correspondingly, we make the replacement e g ^ b ' ^ — » e s ^ b '^F NP (b) m Q 
with the "minimal" ansatz for non-perturbative effects, E112I F NP (b) — exp(— gb 2 ) 
with a non-perturbative parameter g. Fig.l shows the Qt distribution of tDY at 
\/~S = 100 GeV, Q = 10 GeV, y = <f> = 0, and with a model for the transversity 
Sq(x) that saturates the Soffer's inequality at a low scale.EH Solid line shows the 
NLO result using 0, and the dashed and dot-dashed lines show the NLL result 
using ©, F NP (b) = exp(-gb 2 ), with g = 0.5 GeV 2 and g = 0, respectively. 
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